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Abstract 

We investigate the scattering process of impulse in vertical gran- 
ular chain with light impurities. When the perturbation is weak, the 
quantities describing the reflection rate exhibit power law behavior 
with the impuity depth. The exponent is nearly independent of Vi. 
When the perturbation is very strong, the vertical chain shows simi- 
lar behavior to that of the horizontal chain, so the exponent is zero. 
Our numerical investigation begins from the weak perturbation region 
and extend to the nonlinear region and found a peak of the exponent. 
The difficulty in extending the numerical investigation to a stronger 
perturbation region is analyzed. 

PACS numbers: 45.70-n, 43.25. -Fy,46.40.Cd 
Key Words: granular chain, scattering process, power-law be- 
havior, reflection rate 

1 Introduction 

In the last ten years physicists have been interested again in the study of dry 
granular materials. Granular materials are ubiquitous around us. They have 
many properties that are greatly different from those associated with common 
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solids, liquids and gases. These unique properties lead to many applications. 
Many new ideas have emerged from this recent development. These systems 
have been studied from various different viewpoints experimentally, numeri- 
cally, and theoretically [1, 2, 3]. Nesterenko[4]noted that the propagation of 
a perturbation in loaded chain with Hertzian contacts possesses soliton-like 
features. The soliton-like behavior of the signal in two-dimensional granular 
beds was studied by Sinkovits and Sen[5, 6]. The soliton-like properties of 
the signal in the chain of iron balls were examined experimentally by Coste 
et al[7]. Sen et al [8] studied the soliton-like pulse in perturbed and driven 
Hertzian chains and their possible applications in detecting buried impuri- 
ties. For the vertical granular chain under gravity with power-law-typc con- 
tact force, J. Hong et al [9] showed that there are two kinds of propagating 
modes, quasisolitary and oscillatory, depending on the strength of impulse. 
The type of dispersion and the oscillation frequency, wavelength, and period 
follow power laws in depth; J. Hong and Aiguo Xu[10] found that the signal 
speed follows a power law in depth and the exponent varies smoothly from 
a saturant value to zero when the strength of the perturbation varies from 
very weak to very strong. 

To explore the prospects for the application of impulse for detecting cer- 
tain buried objects in granular beds, it is imperative that one must acquire 
a detailed understanding of the its properties. Up to now many problems on 
this topics are still openning. In present paper, we study the scattering pro- 
cess of impulse when it encounters some other kind of grains which we name 
impurities. The granular chain system can be described by the following 
motion equations 



rrinXn = V'(Xn-l, Xn) " V'(Xn, Xn+l) + Ulng, 

n = 1, 2, N. (1) 

where m„ is the mass of the nth grain, V{xn,Xn+i) is the inter-grain inter- 
action between the nth and the(n -|- l)th grains. In this paper we study the 
impulses in a special case with V{xn-i, Xn) — [?7/(p+ 1)] [Aq — (x^ — 
where Aq is the diameter of the spherical grain, p the exponent of the power- 
law type contact force, and rj the elastic constant of the grain under consid- 
eration. Thus equation (1) becomes 
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+171^9, n 



= 1, 2, TV. 



)Y - [Ao - {Xn+l - Xn)]P} 



(2) 



For the Hertzian chain, i.e., the chain with Hertzian interaction between 
neighboring granular spheres, p = 3/2. 



In our system most of the grains are same and the second kind of grains 
only occur as impurities. Hence to analytically get some information of the 
system, we first fix the mass of all the grains to a constant, i.e., m„ = m for 
all n. Then we introduce a new variable ipn, denoting the displacement of 
nth grain from equilibrium, defined by 



where the last term in the right-hand side is the sum of grain overlaps up to 
nth contact and we set Xq — ipo — 0. Eq. (2) is transformed into 



2 Analytic treatment 
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(3) 



1=1 





using Eq. (3). 

For the weakly nonlinear regime, the condition 





is valid and the expansion of Eq. (4) under this condition reads 



'n — 



IJ'ni^n - ^n-l) + /^n+l " ^n) 



(6) 
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where Hn = Tnpg ( j n p is the force constant of nth contact of the 
hnear horizontal chain. We drop weak nonhnear terms of the expansion in 
Eq. (6) to make the system hnear. Thus the vertical granular chain becomes 
a horizontal chain with varying force constants in which the gravity effect is 
contained. Both left and right side of Eq. (6) are linear in V'n- Therefore, 
the scaling analysis tells us that the equation of motion (6) has nothing to do 
with the initial impulse Vi. For this case we[10] have obtained the following 
power-law, 

oc (7) 

The reducing of signal amplitude mainly due to the dispersion effect. 

In our present studies two distances are important. The first one is the 
depth of the impurities dimp- The second one is the distance between the 
impurity center and the grain where we get datas, d. Our main objectives 
are to get some knowledge on the reflection rate and transmission rate of the 
kinetic energy from the impurities. For fixed d, we expect the reflection rate 
follows a power- law and the exponent is independent of Vi. 

For the strongly nonlinear regime, the condition 

i^n-l - V'n) » (8) 

is valid and the expansion under this condition leads Eq. (4) to 

- V [ii'n - i>n+lY + Pgn+l{^n ~ V>„+l)^"^] (9) 

where — {mgn/rif'^^ denotes grain overlap at nth contact. The gravity 
term can be neglected in the highly nonlinear regime, since the gravity effect 
appears in the coefficient gn- Different order of ipn in the left and right 
side of Eq. (9) implies that Wj-dependence must be appeared in the signal 
characteristics. Under the limitation — > oo, Eq. (9) recovers to 

= ?7[(V'n-l - i^nY " {i^n " i^n+lY] (10) 

That is to say, the vertical chain shows same behavior as that of the horizontal 
chain. For the horizontal chain the signal speed is a constant, so the exponent 
of the depth-dependence is 0. 
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3 Simulation results 



Our results are based on careful numerical integration of the coupled equa- 
tions of motion for granular chains with spheres, where varies from 
1000 to 5000 according to our need. The predictor-corrector method, Runge- 
Kutta method, Richardson extrapolation and the Bulirsch-Stoer method[ll] 
are used. The program units are just the same as those in Refs.[8] and [9]. 
That is to say, the units of distance, mass, and time are 10~^m, 2.36 x 10~^kg, 
and 1.0102 x 10~^s, respectively. This kind of program units gives the grav- 
itational acceleration g = 1. Wc set the grain diameter 100 and rj = 5657. 

Figure 1 shows four snapshots for a scattering process, where the impurity 
depth is dimp — 500, which is described by the grain number; there are five 
impurities around the depth dimp — 500; the mass of the impurity grain 
is rriimp — 0.1; the initial perturbation strength is Vi = 0.1. Fig. 1 (a) 
shows a snapshot before the scattering process. Figs. 1 (b) and 1(c) show 
the scattering process. Fig. 1 (d) shows a snapshot where the scattering 
process nearly finished. We use "1st" and "2nd" to denote the leading and 
second leading peaks. The leading and second leading peaks of the incident 
impulse are pointed by an arrow, respectively, in Figs. 1 (a) and 1(b). The 
leading and second leading peaks of the scattered impulse are pointed by 
an arrow, respectively, in Figs. 1 (c) and 1 (d). It is clear that the leading 
peak of the scattering impulse is upwards, which is a characteristic property 
of the scattering process from light impurities [6]. We use five impurities 
because that the scattered impulse is too weak if we use only one impurity. 
To understand the scattering process, we are interested in the variation of 
the ratios vl^j/v}^^^ , v'^^'^/vf^^ and vl''^f / vl^^wiih the depth of the impurities 
dimp-i where v]^^ and u^gj- are maximum values of the leading peaks of the 
incident and reflected impulses, v}^^ and v^^f are maximum values of the 
second leading peaks of the incident and reflected impulses. We obtain these 
datas from a gain with depth n = dimp ~ d. The three ratios describe the 
reflection rate from different sides. During the scattering process the shape 
of the impulse is not stable, so an appropriate distance d is necessary. These 
ratios describe the reflection rate of the kinetic energy in the vertical granular 
chain. In order to get exact value for the leading and second leading peaks, 
we use the smooth velocity-time curve of the measured gain, instead of the 
rough snapshots. 

Fig. 2 shows the variations of vll\/v]^l. , vf^f/vl^^ and v'^^f/vj^^^ with the 
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impurity depth dimp- Fig. 2 (a) is for vlf^ / v}^l{dimp) and vf^f/vf^^^dimp), 
where d = 200, Vi = 1. The squares are numerical results for vll^ vl^J^dimp) 
and the circles are numerical results for vl''^f /vf^^{dimp)- The hues are cor- 
responding fitting curves. Fig.2 (b) is for v^^j / vl^l^dimp) , where d — 200, 
Vi — h. The solid squares are for numerical results and the line is the fit- 
ting curve. We use the simplest power-law model, Vref/vinc oc d~J^p , to fit 
the numerical results. The fitting results suggest that the three ratios follow 
power-laws when dimp is large enough. Prom Refs. [9] we know that when 
the initial perturbation is weak the propagating mode is oscillatory, when the 
initial perturbation is strong the propagating mode is quasisolitary. It is in- 
teresting to check if the power-law model works for both the weak and strong 
perturbations. If it works, we have the following expected result: When Vi 
is large enough, the behavior of the vertical chain is similar to that of the 
horizontal chain; hence there is no power-law. That means /3 — > as — > oo. 

Fig. 3 shows the variation of the exponent for v^^j /v}^^ with the strength 
of initial perturbation, where the line with squares is for the case oi d = 100, 
the line with circles is for the case oi d = 200. The two curves show similar 
behaviors. It is clear that the exponent is independent of Vi when Vi is very 
small. For the case of d = 100, /3 = 0.19 ± 0.01. For the case of d = 200, 
(3 = 0.23 ± 0.01. In fact when is too small we can hardly get helpful 
information from the simulation datas due to numerical errors. The curve 
for d = 100 has a peak around Vi = 5. The curve for d = 200 has a peak 
around Vi — 10. When Vi > v , (3 decreases quickly with the increasing of Vi. 
We tried to add more points to the figure, but when Vi is larger than 20 we 
can hardly measure the second leading peak from the incident impulse. We 
can know the reason from Fig. 4. 

Fig. 4 shows a f„(t) relation, where dimp = 300, d = 200, n = 100, 
Vi = 15, and rriimp = 0.1. The impulse near t — 2 is the incident one. The 
impulse within 9 < i < 10 is the reflected one (Only part of the reflected 
impulse is plotted). It is clear that for a strong perturbation the propagating 
impulse is quasisolitary and it has no evident the second leading peak; at the 
same time, the leading peak of the reflected impulse is very small and there 
exist some oscillations before the leading peak occurs. 

When the perturbation is weak, the propagation of the wave is oscillatory; 
hence the two ratios v^f^/vl^l and v^^f/'^^inc describe well the reflection 
rate of the kinetic energy from the impurities. But with the increasing of 
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t'i, the propagation of the wave is more and more quasisohtary and the ratio 
'^refhlnc describes the reflection rate better. Fig. 5 shows the variation of 
the exponent (3' for vl^j / vl^J^dimp) with Vi, where d — 200. It shows similar 
behavior to the exponent (3 for v^^f /vf'!^^{dimp)- When Vi is very small, (3' w 
0.2. 

Although we have found the power-law behavior for v^fj / vl^^^{dimp) for 
not very large Vi (see Fig. 2 for an example), it is hard for us to get a figure 
which is similar to Fig. 3 and Fig. 5. The reasons will be given in the following 
part. 

For a fixed distance d, when the depth of impurities dimp is large, we 
find the curve vlf^ / v}^\{dimp) oscillates greatly with di^p. Fig. 6 shows two 
examples, where d = 200; = 6 in (a) and = 4 in (b). We can understand 
the reason for this kind of behaviors through Fig. 7. Fig. 7 shows several 
Vn{t) curves from which we can study the scattering process of the impulse, 
where n = dimp - 100, = 0.1, rriimp = 0.1; di^p = 600 in (a), dimp = 800 
in (b), dimp = 1000 in (c), and dimp = 1100 in (d), respectively. For the 
refiected impulse, only the leading peak and the second peak are shown. For 
a fixed value of d, when the depth dimp is small, the refiected impulse and 
the incident impulse are separated, i.e. there is a time during which the 
measured grain is static. For these cases, we can get exact values for vj^ll, 
v^gj, vf^^, and f^ej. When the depth dimp is large enough, the tail of the 
incident impulse and the refiected impulse are connected. For these cases, we 
can hardly get exact values for v^f^ and f ^g^. From Fig. 7 we also understand 
the fact these phenomena first infiuence the measurement of v^fj , then v^"^. 

We name the increasing process of dimp before a power-law occurs a tran- 
sient process. The simulation results show that the ratio v^^} / v}f^l{dimp) needs 
a much longer transient process than v'^.'^j / v'j^^{dimp) ■ The initial impulse is 
stronger, we need a longer transient process. The curves for v'^^f / vf^^ {dimp) 
decreases with the increasing of dimp, while the curves for vll^ v}^l{dimp) first 
increases, then decreases. The initial impulse is stronger, the curves decrease 
more slowly. We also notice that there are more fitting errors when Vi is 
large. 
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4 Conclusion 



In present paper we study the scattering process of impulse in vertical gran- 
ular chain with light impurities. Most of the studies resort on numerical 
simulations and the analytic treatment helps to understand the behaviors. 
When the perturbation is weak, the quantities describing the reflection rate 
exhibit power law behavior with the impuity depth. The exponent is in- 
dependent of Vi- When the perturbation is very strong, the vertical chain 
shows similar behavior to that of the horizontal chain, so the exponent is 
zero. Our numerical investigation begins from the weak perturbation region 
and extend to the nonlinear region and found a peak of the exponent. The 
difficulty in extending the numerical investigation to a stronger perturbation 
region is analyzed. More work needs to be done to completely understand 
the scattering process of the signal in vertical granular chain. 
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Figure caption: 

Figure 1: Snapshots for a scattering process, where the impurity depth is 
dimp — 500, which is described by the grain number; there are five impurities 
around the depth dimp = 500; the mass of the impurity grain is rriimp — 0.1; 
the initial perturbation strength is Vi = 0.1. 

Figure 2: Variations of vll^vj^^ , v^^f /vf^^ and v'^l'^/vj^l with the depth 
of the impurities di^p- Fig.2 (a) is for vll) / v}^l{dimp) and vl'^f /vf^^{dirap) 
where d — 200, Vi — 1. The squares are numerical results for vliy vl^J^dimp) 
and the circles are numerical results for v^^f /vl^cidimp)- The fines are corre- 
sponding fitting curves. Fig.2 (b) is for v'^'^f / '^inlidimp) where d = 200, Vi — 5. 
The solid squares are for numerical results and the line is the fitting curve. 

Figure 3: Variation of the exponent for vf,l''j / vfj^^ with the strength of 
initial perturbation, where the line with squares is for the case of ci = 100, 
the line with circle is for the case of d — 200. 

Figure 4: The Vn(t) relation, where di^p — 300, d = 200, n — 100, 
Vi — 15, and mimp = 0.1. The impulse near t = 2 is the incident one. The 
impulse within 9 < t < 10 is the refiected one (Only part of the refiected 
impulse is plotted). The second leading peak of the incident impulse and the 
leading peak of the refiected impulse are very small. Some oscillations occur 
behind the main incident impulse. 

Figure 5: Variation of the exponent (3' for v^'^f / v}^\{dimp) with Wj, where 
d = 200. It shows similar behavior to the exponent P for v^2f /''^incidimp)- 
When Vi is very small, /?' ~ 0.2. 

Figure 6: Refiection rate vl^j/vj^l versus di^p, where d — 200; Uj = 6 in 
(a) and Vi — A in (b). 

Figure 7: Velocity Vn{t) curves, where n = dimp — 100, Vi = 0.1, rriimp = 
0.1; dimp = 600 in (a), dimp — 800 in (b), dimp — 1000 in (c), and dimp = 1100 
in (d), respectively. 
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